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Abstract 

The contribution of scale-dependent vacuum expectation values (VEVs) of Higgs 
scalars, which gives significant effects in the evolution of fundamental fermion masses 



Oh! 

r-| ! in the Minimal Supersymmetric Standard Model (MSSM), is now considered in the 



derivation of the one-loop analytic expression for the evolution of the left-handed 
Majorana neutrino masses with energies. The inclusion of such effect of the running 
VEV increases the stability of the neutrino masses under quantum corrections even 
for low values of tan/5 > 1.42 at the scale /x = lO^^GeV, and leads to a mild decrease 
of neutrino masses with higher energies. Such trend is common to that of other 
fundamental fermion masses. 
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In recent years a large number of theoretical papers were devoted to building 
models for generating small neutrino masses and lepton mixings within or outside the 
framework of the Grand Unified Theories (GUTs) with extended U{1) group[l]. Both 
analytic and numerical studies [2, 3, 4] have been carried out for checking the stability 
of the textures of neutrino mass matrix and lepton mixing matrix under quantum ra- 
diative corrections [5]. There are basically two approaches: the top-down approach [2] 
which predicts the neutrino masses and mixings in terms of GUT-parameters, and 
the bottom-up approach [6] which predicts the running parameters at higher scales 
in terms of experimentally determined values at low energies. In the top-down pro- 
gramme, one usually starts with the running of a set of the RGEs for Yukawa matrices 
and gauge couplings in the MSSM (or SM), with three right-handed heavy neutri- 
nos, taking into account the effects of the heavy neutrino mass thresholds, from GUT 
scale down to the lightest right-handed neutrino mass scale {Mf>i). This fixes the left- 
handed Majorana neutrino mass matrix ^^^(Mri) through the see-saw mechanism[7] 
at this scale, 

rriLLiMm) = vlY,{MRi)M^},Yj{Mm) (1) 

Below this scale Mjn the right-handed neutrinos decouple from the theory, and the 
neutrino mass matrix in Eq.(l) is taken as [2] 

mLL(Mm) = vIk{Mr,) (2) 

where k, is the coefficient of the dimension 5 neutrino mass operator. In the energy 
range from Mjn down to low energy at rrit, the running of the coefficient k in the 
diagonal charged lepton basis, fixes the neutrino mass matrix at scale rrit, 

rriLLirnt) = vlK{mt). (3) 

In the above discussion only the scale-dependence of k is considered, and not the 
running of the vacuum expectation value(VEV),fn in Eqs.(l)-(3). This led to the in- 
crease of neutrino mass eigenvalues with energy scales, giving significant effect for low 



tan j3 values. As it is strongly tan (3 dependent, this efect may lead to the instability of 
the neutrino masses under quantum radiative correction^. For higher values of tan (3 
the stability is again improved. Such increasing trend of neutrino mass eigenvalues 
with the increase in energies, is opposite to that of the general trend shown by other 
fundamental fermions (charged leptons and quarks) [8,9]. The effects of the contribu- 
tions of the scale-dependent vacuum expectation values(VEVs) of Higgs scalars in the 
one-loop analytic expressions, in the evolution of quarks and charged leptons masses 
at higher energies in the MSSM, had been studied in Ref.(8), and this effect is quite 
significant. 

In this paper we study the stability of the magnitudes of neutrino masses at low 
tan/3 and their running behaviour at different energies, by considering the scale- 
dependence [10] of the vacuum expectation value (VEV), Vu, along with that of k. 
The expression in Eq.(3) is now modified as 

rriLLit) = vlityit) (4) 

where Vuito) = fosin/3, vq = 174GeV, t = In/i, to = Inm^. The above equation(4) 

can be written as 

dlnmLLJt) _ dhiK{t) dlnVujt) 

dt ~ dt dt ^ ^ 

where the second term on the right-hand side of the above equation is the contribution 
from the running of the VEV. The RGEs for Vu [8,10] and k, [2,5]in the diagonal 
charged lepton basis, for one-loop order in MSSM, in the energy range t > to, are 
given by 

[7^^1 + 7^2-3/ii], (6) 



and. 



dt 167r2'20 ' 4 

dlnK 1 r6 



dt 167r2^5 



[-gf + Qgi - Qhf - 5,,K ' h,h'^ (7) 



^In Refs.[3,4] the stability condition is decided by the change in texture of neutrino mass matrix 
only. Here we emphasise that changing pattern of the overall magnitudes of neutrino mass eigenvalues 
at different energies, may also cause instability. 



respectively. Substituition of Eqs.(6,7) in Eq.(5) gives 

Upon integration from low scale to = Inmj to high scale tm = InM^ji where tjn > to, 
we get the correct expression for the neutrino mass matrix at to, 

{jnLL{tRi))ij 

1 /•lnJ\/m 



167r^ Jlnmt \ gi{to) J 

where f = t,T; i = 1,2,3, and bi = (33/5, 1, —3) for MSSM. The correct expression 
in Eq.(9) will certainly affect the earlier numerical results obtained without taking 
the effect of the running VEV [2] at scale M^i . 

For simplicity we now follow the analysis of the RGEs for neutrino mass eigen- 
values[4]. With the inclusion of such scale-dependence of VEV in Eq.(6), the RGEs 
for the mass eigenvalues given in Ref. [4] in the diagonal charged lepton basis, is now 
modified as, 

amnii^a t .^-^ r 9 2 ^ 2 , o!,2ta2i /ioA 

-^r- = ^\^ f-i5^^ - 2^^ + '^^^^«] ^^'^ 

where a = 1, 2, 3, and Vf,a is the MNS mixing matrix element. The correct expression 
for the neutrino mass ratio at different energy scales, is also obtained by integrating 
Eq.(12) as 

In getting Eq.(13) we have neglected very small effects due to J^^e compared to /^, 
and also assumed K-a does not change much in the integration rangeQ. For a typical 

value of the element of MNS mixing matrix Vrz — "To' ^^ *"^^ ^st the condition 

■^ Such approximation can be justified for the calculation of the mass eigenvalues and their ratios 
as the second exponential term in Eq.(13) gives almost 1 for low values of Ir- 



iTT-usito) > iTT-usitRi) following Eq.(13), which shows a mild increase in neutrino masses 
with the decrease in energies, even for small tan/? > 1.42. This is due to the fact that 
the ratio Rsitm) is now independent of e^'^* in the first exponential factor in Eq.(13). 
The same is true in Eq.(9). In fact the contribution of the running VEV effectively 
makes the replacement in the exponential factor: 

in Eqs.(9,13). 

We now study the effect of the running VEV in the evolution of squared neutrino 
mass difference, A-mf- = |m^j — ""^Ll with energies. By taking square on both sides 
of Eq.(13), and considering for two mass eigenvalues a = i,j, we get approximately, 

A m^^.(tm) ^ AmJ,.(to)e*='(^^.i+f^«2)e4^"^- (15) 

where we assume that the small difference between Vri and Vrj for i,j = 1,2,3, 
does not alter much the last exponential term which can be approximately taken as 
e^^ri^T ~ e^ ^i ^ ^ 1 for the low values of J,-. This amounts to neglecting small changes 
in the texture of neutrino mass matrix which would be relevant for the evolution of 
mixing angles. The evolution of Am'fAtjii) is now stable with the effects of running 
VEV for both low and high values of tan f3, otherwise it would have been more strongly 
tan/9-dependent with e^^^* in the exponential factor in the case where the effect of 
running VEV is not included, causing more instability at low tan/3 values. 

The running of the ratio of two neutrino mass eigenvalues, R23 = Tny2/fnu^ (and 
hence the running of RR23) is independent of the effect of running VEV, so that the 
ratio of ratios is, 

RR2s{tn.) = ^^ - e-^^-^^ (16) 

where 

6V,l2 = K'3 - Kl (17) 



which can be either positive, negative or zero. For the positive value, 5V^.^2 > as in 
hierarchical case [2], one gets the condition, 

i?23(to) > R2z{tRl) (18) 

which implies the increase in the neutrino mass ratio 171,^2 /itt-us with the decrease in 
energies[2]. If we start with degenerate neutrinos, 171^2 = f^uz at the scale M^i, then 
we would get mi^2 > f^v^ at scale m^. This shows that nearly degenerate neutrinos 
are not stable under quantum corrections [3]. 

The above relations in Eqs. (16-18) for a = 2, 3, can be generalised for any pair of 
mass eigenvalues a = i,j. For inverted hierarchical case [2] with mj,i > 171,^2, we may 
have (5^4^21 < which leads to 

Ruito) < Rl2{tRl) (19) 

where the neutrino mass ratio 771,^1/771^2 decrease with the decrease in energies[2]. The 
effect of the running VEV does not change the textures of neutrino mass matrix and 
hence the MNS mixing matrix. 

Next we turn to numerical analysis of the RGBs in the bottom-up approach in 
running from low energy scale to to high energy scale, replacing tm by running t in 
the above equations (9-19). We make use of the following input values of the running 
fermion masses 7ni{7ni) of the third family: 

mt,b,r = (166.5, 4.2, l.785)GeV (20) 

where, for heavy flavours(top and bottom quarks) the values are derived from input 
pole-masses 712^°'^ = 175.6GeV [11] and 7nl"'^ = 4.7GeV [12,13] using two-loop RGBs 
in QGD. The initial input values for the top, the bottom and r-lepton Yukawa cou- 
plings at top-quark mass scale to = Inmj in the RGBs in MSSM, are usually obtained 
as 

ht{to) = 7nt/{174smp),h,r{to) = 7nb,r/{l747]h,rCosp). (21) 



Using the CERN-LEP measurements at Mz = 91.18GeV, 

a3{Mz) = 0.118±0.004,a-^(M^) = 127.9±0.1, sin^ ^^(M^) = 0.2313±0.0003, (22) 

we obtain the values of gauge couphngs at scale to using one-loop RGEs, assuming 
the existence of one-light Higgs doublet (n=l) and five quark flavours below rrit scale, 

a^l^ito) = 58.42, 29.67, 8.89 (23) 

The QCD-QED rescaling factors [6] are calculated as 

77/ = (1.54, 1.017),/ = 6, r (24) 

As a result of the numerical analysis of the RGEs for Yukawa and gauge couplings 
at two-loop level [6] in the energy range tg < t < tjj, the unification of three gauge 
couplings is observed at Mjj = 1.82 x lO^^GeV. The values of Yukawa couplings 
{ht,hb,hr), gauge couplings and values of integrals Jj defined in Eqs.(10,ll) for dif- 
ferent values of tan/? = 1.42 — 60.0 are estimated at different energy scales. 

We present our numerical results in Figs. (1-4) where the solid line refers to the 
analysis with the effects of running VEV in the present calculation (referred to as 
case A). We also present the corresponding results without the effect of running VEV 
in dotted line (referred to as case B) for comparison only. With a typical input value 
Vr3 = 4^, the variation of the ratio Rsit) defined in Eq.(13), with energy scales t for 
two representative values of tan/5 = 1.63 and 57.29, are presented in Figs.l and 2 
respectively. These figures show the evolutions of neutrino mass eigenvalue m„3 with 
the increase in energy scale. 

We observe that for high value of tan/? = 57.29 in Fig.2, the evolution of the ratio 
i?3(t) = 772,^3 (t)/m,^3 (to) is almost stable in both cases A and B. However, for low value 
of tan/3 = 1.63 in Fig.l, there is a significant increase in -R3(t) at higher energies in 
case B. For example, at /i = 1.82 x 10^^ GeV, the ratio -R3(t) is about 5.63 in case B 



as shown in Fig.l by dotted line. Such unwanted feature which may cause instabihty, 
is not present in case A (sohd hne in Fig.l). Fig. 3 shows the variation of neutrino 
mass, Rsitjii) = m^^itjii) / my^ito) at a particular scale, tjn = 27.63 corresponding 
to Mp>i = 10^'^GeV, with different values of tan/3 = 1.42 — 60. We see that at low 
tan/3 > 1.42 region there is a significant enhancement in R^ < 5.3 in case B whereas 
the ratio is stable in case A for all values of tan/3. For higher values of tan/? the ratio 
is again stable in case B. The same analysis is true for the cases of other two mass 
eigenvalues with a = 1, 2. Similar analysis can be done for the evolution of A-m?- in 
Eq.(15), which would be very unstable in the low tan/3 region in case B. However, it 
is now stable for all values of tan j3 under radiative corrections at higher energies in 
case A. 




Fig.l Variation of R^it) = 'Ti^3(t)/?Tiiy3(to) with energies t = In/i for small value of 
tan/3 = 1.63. The results with and without the effect of running VEV, are shown 
in solid line and dotted line respectively. 




Fig. 2 Variation of i?3(t) — muzify j rtivziti^) with energies t = In/i for large value of 
tan/3 = 57.29. The results with and without the effect of running VEV are shown 
in solid line and dotted line respectively. 
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Fig. 3 Variation of i?3(ii?i) = m^3{tRi)/m^3{to) with ln(tan/3) for Mri = 10^'^GeV. 
The results with and without the effect of running VEV are shown in solid line and 
dotted line respectively. 
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Fig. 4 Variation of the ratio of the neutrino mass ratio RR23,{t) — (m(X )/{ "m(t\ ) 
with energies t = In/i for large value of tan/3 = 57.29. 

Finally, we study the relative rates of the evolution of two neutrino mass eigenval- 
ues in terms of their ratio, i?23 = 'nT'u2 / itt-us given in Eq.(16), in going from low to high 
energies. We consider high value of tan/5 = 57.29 where the effect of /,- is large, and 
this ratio increases with the decrease in energies by a few percent only. This is shown 
in Fig.4 where we present the evolution of the ratio of the ratios RR2z{t) in Eq.(16) 
with energies. This leads to a mild increase in the hierarchical relation, 171^2 /it^-us at 
lower energies. As noted earlier, such hierarchical ratios are independent of the effect 
of running VEV. Finally we point out the changes aring from the running of VEV 
in the earlier calculations [2] of neutrino masses. The earlier results at low scale rrit 
in Ref.[2] do not change at all. However if we prefer to express the neutrino mass 
matrix at higher scale M^i, then we have to take the effect of running VEV, f„(t/ji) 
in place of ^^(to), which modifies the earlier numerical results at the scale Mm. 



To conclude, we have considered the contributions of scale-dependent vacuum 



expectation values (VEVs) of Higgs scalars in deriving one-loop analytic expression 
for running the left-handed Majorana neutrino masses with energies in the MSSM. 
This gives significant changes in the expression of the evolution of neutrino masses, 
and also increases the stability of the neutrino masses under quantum corrections 
even for low tan/3. We observed a mild decreasing trend of neutrino masses with 
higher energies, which is now common to that of all other fermion masses in nature. 
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